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Light quark confinement|^ Q 



We have described the confinement of heavy quarks in an analogy with the theory 
of the supercharged nucleus [1,2]. Let us now suppose again that a is behaving 
like 

a 




Fig. 1. 

Making this assumption, we are neglecting gluon-gluon interactions and the ex- 
istence of gluons as real particles. Our aim is to see, what can arise from the 
discussion of light quarks only. We introduce A corresponding to Uc and consider 
quark masses mg <C A. The interactions of light quarks (for which ttiq <C A) will 
be discussed in a rather simplified way. We will take into account all possible 
interactions 



where the gluon propagator (considered as an effective photon), corresponding to 
the dotted line is 



D 



a 



Further, we look for a model which enables us to see, what happens to the 
fermions if there is an interaction between them, as indicated above. The question 
is, how the bound states or the Green function behave in such a case. 

Let us consider the energy of two quarks, u and d, for example. Without 
interaction there will be positive energy states with E > 2m and negative energy 



^This is the third lecture on quark confinement given by V.N.Gribov in 1992 in Orsay. An 
extensive discussion of the consequences of all this for the structure of the Green function can 
be found in [5,6] - in the two last papers concluding his 20 years long study of the problem of 
quark confinement in QCD. 

^The text was prepared for publication by Yu. Dokshitzer, B. Metsch and J. Nyiri on the 
basis of a tape recording and notes taken during the lecture 
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states with E < —2m: 

E{ud) 



2m 



a 



-2m 



Introducing the interaction, for smaU a we will find that there are some bound 
states near 2m and —2m. 
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So far, we consider the usual Dirac vacuum: the negative energy states are oc- 
cupied, and the positive ones are empty. Increasing the coupling, i.e. increasing 
a, we could expect that the magnitude of the energy is decreasing and the levels 
corresponding to the bound states will come closer and closer to zero. 




-2m 



a 



With a further increase of the coupling up to a critical value, one possibility for 
the levels will be just to approach the zero line and never cross. There is, however, 
also a possibihty of crossing. We will see that the first case corresponds to normal 
spontaneous symmetry breaking. But, if the levels cross, and especially in the 
most clear case, when they pass the lines 2m and —2m, respectively, everything 
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will change and we arrive at very different phenomena: 
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Indeed, now the original vacuum which corresponds to the case when level 2 is 
empty and level 1 is occupied, is absolutely unstable. We have to fill the new 
negative energy state and leave the positive energy level empty. But by filling 
this new state, we get an excitation, a meson-type state with a mass /i. For 
free quarks this would mean that the quark with negative energy decays into a 
negative energy meson (filling the negative energy levels) and creates a positive 
energy quark. As a result, the Dirac picture in which all negative energy levels 
are filled up and all positive energy levels are empty, is destroyed. But if so, a 
positive energy quark also decays into a positive energy meson and a quark with 
negative energy. This means that both decays 



are possible, and both q~ and g+ are unstable. 

The question is now, how to deal with the bound state problem. Of course, we 
could just start to calculate the bound states, considering the interactions without 
corrections to the Green function. However, one has to take into account that 
the fermion-fermion interaction changes the effective mass of the quarks and this 
in its turn will change the bound states considerably. 



which makes the problem more complicated. We thus will have to consider bound 
states and the Green functions on equal footing. 

Up to now, the only approach which deals with this problem and is self- 
consistent is the Nambu - Jona-Lasinio model [3] . It considers the fermion Green 
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function corrections due to a four-fermion interaction: 




In spite of the strong dependence on the cut-off, the model preserves all sym- 
metries in the Green function and in two-particle interactions. Let us present 
the result of Nambu and Jona-Lasinio in a way somewhat different from what 
is given in [3]. We express it as the dependence of the renormalized mass m on 
the bare mass ttiq. They found that if the effective coupling (it depends on the 

2 

definition in their case) |^ is less than unity, the curve will be just the usual one: 



m 




If, however, the coupling ^ is larger than unity, the dependence will be like Q: 



m 




Fig. 2. 

•^This result is not always quoted, but it is present in their paper. 
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According to the interpretation of Nambu and Jona-Lasinio, the upper part of 
the curve, which at rriQ = reaches a finite point, corresponds to the spontaneous 
symmetry breaking. But there are three solutions at ttlq = and at sufficiently 
small mo values. What Nambu and Jona-Lasinio claim is that the lower part 
of the curve is unstable, and there is a real vacuum. I agree with this, if mo > 
TTT'crit- Wc Can ask: what is the source of instability of this curve? The general 
argumentation is the following. Talking about spontaneous symmetry breaking, 
m is like a magnetic field in a ferromagnetic; we just choose a definite direction. 
But mo is like an external field, and the system is like a compass. If the external 
field and the induced field are pointing at the same direction, the situation is 
stable. If they point to opposite directions, the compass will change. 

I am, however, not sure that the instability of the almost perturbative solu- 
tion which contains no condensate at all can be explained in such a way. The 
explanation may be right for the part b of the curve in Fig. 2 which corresponds 
to a big spontaneous magnetic moment and the opposite direction of mo. It does 
not work for the part a close to perturbation theory which has no spontaneous 
magnetic moment. And, looking more carefully at the curve, we see that the 
part b corresponds to pseudoscalar states inside the Dirac see, while on the piece 
a both the pseudoscalar and scalar states are inside, both levels passed. Recog- 
nizing this, one can conclude that indeed, the mentioned state is unstable, but 
for a trivial reason: the corresponding level is inside the Dirac see and it is not 
filled up. The problem is, what happens if we fill Tip this level. It remains an 
open question, what can be considered as a ground state under these conditions. 
And it is a problem how to get these results in a more self- consistent way, not 
depending on the cut-off so strongly. 

The Nambu - Jona-Lasinio model can be reproduced in our picture. For this 
purpose, just as a theoretical exercise, let us use a not going to unity, i.e. draw 
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Fig.3. 

instead of the curve in Fig.l. In this case there will be a second scale A2, and 
outside this scale we will have just point-like interaction. This reminds of the 
Nambu - Jona-Lasinio picture, which, apparently, can be reached somehow in 
our approach. The problem is, how to write constructively the corresponding 
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equation, and whether this can be done at all. Of course, this constructive part 
can be only approximate. But if we recognize that it can be written, then we 
will be able to develop a theory in which we put the main ingredient of our 
discussion as an input into our solution and try to find the real construction. 
The main difference will appear in the analytic properties of the Green function. 
The Green function of a fermion for such a case would be quite different in its 
analytic structure compared to the usual one. 

I am afraid 1 will not have the time to come to this point today, but I would 
like to explain just the physics. 

How to write the equation? What happens in the real case and how to deal 
with it? Let us start with with the Green function. What do we know, what are 
we supposed to know about this Green function as a function of (f" ? 



Beyond a certain A in the region where the coupling is small, it is asymptotically 

free; here the Green function has to satisfy the renorm group equation. But if 
as a result of the interaction a mass is acquired, this mass would be somewhere 
at smaller g^; here the equation becomes essentially very complicated and we are 
not able to extract a reasonable structure. 

The idea to write an equation which is correct in both regions, near the 
threshold and at large g^, and to match these two solutions, comes from the 
following consideration. Suppose that acrit/'^ is small: 



Now, however, we may ask: how could new masses, new solutions etc. appear 
at all at such a small a. Obviously, this a has to be multiplied by something 
large. What happens, for example, at large ql We know, that there is always a 
logarithm of g^/A^ and the real parameter becomes 



which is, in spite of the smallness of a, big enough to change the Green function 
essentially. But near the threshold there is also a logarithm: 



m 
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q — m 



with some scale A = A or m 



A2 
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which is always present. In other words, in this region there could be also a 
quantity which changes seriously in spite of the relative smallness of a. Hence, 
we want to write the equation which is correct near the threshold, taking into 
account correctly the singularity of a supposed mass, and after that compare this 
with the renorm group equation; we shall see whether it is possible to write an 
equation which is correct in both regions, and if yes, we will try to solve it. In 
order to get the singularity correctly, we take the second derivative of with 
respect to the momenta. 



(g-</'P__ 

/ \ / / \ \ \^) 
r< — \ 1 I ( S _L -d L \ i_ 

^ -^^^ q G{q') ^ 1' Q" 

The contribution of the first term is trivial, the second derivative of g + m gives 
zero. Taking the second derivative of the first graph, it can be easily seen that 

d'- = -AniS^'\q - q'). (3) 

(g — q'Y — le 

This gives for the first diagram ^^G{q)^n^ - just by direct calculation. In other 
words, we make it local. From this diagram we take the contribution where 
k = {q — q') - the momentum of the photon - essentially equals zero. 

Let us now look at the second diagram. We have here the choice of taking 
the second derivative at one of the photon lines, or to differentiate once at one 
line and once at the second line. Having in mind that all the integrations would 
have a structure which need some logarithmic enhancement, it would mean that 
the most important regions of integration in this integral would be those where 
ki <C k2 and k2 <C ki. We take the derivative at ki and then set it to be zero, 
but for k2 the integration will give the same as before. If this integration gives us 
two logarithms, we kill one and recover it after the integration of our differential 
equation; but we still have the first one. But if we differentiate once one line and 
once the other, we will always sit on the region ki ~ k2, because they have to 
be of the same order. And, in this case, there is no logarithm at all; after the 
integration, we will recover one, but one order will be lost. Clearly, a possible 
approach is to try not to choose different diagrams, but to use the small k region 
of integration. Ordering the integration inside the diagram in such a way that one 
momentum is much smaller than the others, and differentiating this diagram, we 
will find a relatively simple answer. Indeed, suppose that we have any diagrams 
with any loops. If we differentiate some lines twice (it can be any line) and neglect 
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all first derivatives, we get an amplitude of the following structure: 



k=0 k=Q 




This is just the Compton scattering of a zero momentum photon k = 0, and for 
this quantity the most singular contribution is obviously 

r^(0,g)G(g)r^(0,g) 

which corresponds to the diagram 

G{q) 





But the vertex F is at zero momentum and hence r^(0,g) = d^G ^(g). In 
this approximation we can write a very simple equation: 

d'G-\q) = ^d,G-\q)Giq)d,G-\q) (4) 

TT 

which differs essentially from any Bethe-Salpeter type equation. Indeed, using a 
Bethe-Salpeter type equation, we do not change the vertex part and end up with 
rather bad properties. Equation is scale invariant, it is 75-invariant, it has 
many nice symmetry properties and, what is most important, it has a correct 
behaviour near the threshold. 

The gauge is fixed, because we used 

D = ^ 

It is an important question, what we would get in different gauges. In Feynman 
gauge we are very lucky: we find an expression which does not depend explicitly 
on the expression for the Green function. Using a different gauge, we would find 
the infrared behaviour of this diagram to be more complicated, and we would 
not be able to extract universally the region of small k. We would have integrals 
over q which are also possible to use, but with the necessity to think about the 
behaviour near the threshold. 

We, however, have chosen this gauge; we did not destroy the general features 
and used the current conservation which just corresponds to = d^G~^. Ac- 
cepting this, we can now ask, what is the relation to the renorm group equation. 
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Suppose that we would like to write the renorm group equation in the same 
spirit. Let us take again the second derivatives. In this case we would be definitely 
correct, because we know that it is a logarithmic approximation. 

In our logarithmic approximation we will do exactly the same with the only 
difference that a would be a{q^) . In the renorm group equation a is a function 
of q^. But, of course, a in general depends on two momenta: k"^ and (f. And in 
the renorm group equation at large momenta, in the ultraviolet region, a depends 
on the variable which is the largest. Since we are close to A; = 0, this means that 
here we will have ot{(f), and we will recover the renorm group equation at large 
(f. If we solve this equation with a slowly varying a, we will be correct in both 
the threshold region and the ultraviolet region. 

We also have to formulate an equation for the bound states under the same 
assumption. Looking for bound states, we consider scalar and pseudoscalar ver- 
tices. This vertex 




has to be equal to 




+ 





+ 



Here (p{q,p) depends on p, the total momentum of a pair, and q, the quark 
momenta being given by q + p/2 and q — p/2. With the same procedure as in 
obtaining the equation for the Green function we find for the vertex (see [4] for 
some details): 



a 



d (p{q,p) = - A^{q)d^^{q,p) + d^^A^{q) - A^{q)^{q,p)A^{q) 

TT L 



(5) 



where A^ — d^G~^G, A^ — Gd^G'^ . It means that we have two equations in 
this approximation. We used this approximation just to be constructive and to 
study what will result if we make this approximation. In principle, solving both 
equations we will get everything what is necessary: we know G and A^, we have 
a linear equation for bound states, we can see what is the type of the energy etc. 
The equation for the bound states has very nice features. It is beautiful from the 
point of view of the Goldstone theorem in the following sense. 

Suppose I have some symmetry in my equation, e.g. 75-invariance. Since there 
is no 75 in equation (5), it is 75-invariant. But of course the boundary condition 
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for at g — > oo is just Gq^ — (q — nT^o) , ^-nd thus destroys the symmetry. But 
suppose that itlq = 0. In this case there would be symmetry here, which means 
that the Green function will not be unique, since it can be 

G-'^ + 5G-\ 

where G~^ is some solution and SG~^ oc 75^"^ . This means that the variation 
6G~^ also is important. What would be the equation for the variation? If we 
calculate the variation on both sides of equation (4) we obtain 

{SG-') = - (d^{6G-')Gd^G-' + d^G-'Gd^{6G-') - d^G-'G{6G-')Gd^G-') , 

SO we find that 9? = bG^^ fulfils equation (5) at p = 0. It means that if some 
symmetry is broken, i.e. if there are multiple solutions of the equation for the 
Green function, we always will have some solution of the equation for the vertex 
at p = 0, which is the Goldstone. 

It is clear, that in the present model we can discuss many questions, use a 
running coupling a as in Fig. 1 and reproduce the NJL-features without any 
essentials depending on a cutoff. Before discussing this point further, we will 
first look for the solution of (4) and discuss the result. 

Above, we introduced = {d^G^^)G^ which is a very useful quantity. Since 
G~^ = a| + 6 is essentially a 2 x 2 matrix, is a [/(2) gauge potential: 

d'^G-' = d,{{d,G-^) G G-') = {d^A^)G-' + A,{d^G-') = -A,d,G-' (6) 

TT 

where in the last step we used Eq. (4). Multiplying from the right by G we thus 
find 

(X 

9iiAa + A^Aa — AuAf^. (7) 

TT 

This means that 

df^A^ — —(5A^A^] 

and A^ is a pure [/(2)-gauge potential with a condition (3 — 1 — ^. Of course, 
this is just a useful trick. Important is to write down the real equation for the 
Green function. The most natural thing is to express in the form 

where n is a 2 x 2-matrix 




It is just easier to use this form for our purpose: we can find an equation for p 
and an equation for ip. Both are functions of q^: p{q^)i ^{q^)- There is, however, 
no scale in the equation; it contains only a derivative of q. If we introduce 
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we will find an equation in which ^ can be considered as "time", and which is an 
oscillator equation. In fact there are two oscillators, one for p, the other for (/?, 
and they will satisfy non-linear equations. For (/? we find 

<^ + 2 + 0-3sinh(^ = O. (8) 

This is just an oscillator with damping; a similar equation can be written for p. 
Important is that that there has to be "energy" conservation in this equation. 
Indeed, we said that ^ plays the role of time; it, however, did not enter the 
equation explicitly. Thus there has to be a conservation law which, as it is easy 
to show, leads to 

(l+/3g' = l+/3^(^-3.inh=f). (9) 

We thus can eliminate p altogether, and find the equation for tp 

(p + 2, 



-/32 |^3sinh2|-^j0-3sinh99 = O^ 



which is an oscillator with damping. Having this in mind is sufficient to under- 
stand the structure of the solution. Indeed, what is this </? ? We have 

G-^^pcosh^ + ^psinh^. (10) 

The perturbative solution is if close to iir. In this case the first term is zero, the 
other is proportional to q/q - this corresponds to the massless solution. Since mo 
is small, we have to have solutions like this at g — oo. 

Now we have to find the solution everywhere. Let us first investigate the 
equation without damping; we get 

(/3 — Ssinh^? = 0. 

If we go to the Euclidean space, (f — iip, the potential becomes a periodical 
potential: 

ip — Ssinip = 0. 
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7[ ^TT 




Fig. 4. 

We have to look for a possible solution for this structure with damping. What 
does this mean? For the oscillator with damping any solution at ^ — > oo has to 
be in a minimum, because the energy is decreasing. But if ^ is going to — oo, the 
energy is growing. What could be in this case a normal, reasonable solution? It 
is almost clear that the only possibility is to put at ^ — > — oo the "particles" in 
this oscillator at the maximum, and start to move them slowly; eventually, they 
will appear inside the well. 

There is a most important question, namely: what is the critical coupling in 
this case? What do we know about an oscillator with damping? If the damping 
is large enough, all the trajectories will go monotonically to the minimum. If 
the damping is sufficiently small, the solution will start to oscillate. In order to 
see when this transition happens, we have to look for the equation just near the 
minimum ijj — n. With = tt — ■0 we have for small (p 

4> + 2^1 + 3(3''<j) + 30 = , 
with fundamental solutions 

0i_2 = e"'''^ where 1/1,2 = —/iTs^ ± ^3/32-2. 
So we have monotonic behaviour for 3/3^ — 2 > 0. On the other hand if /S^ < |, 

we will have oscillations before reaching the minimum. The critical angle ipc, 
which separates the regions where the solution is monotonic and where it oscillates 
can be shown (see e.g. (4)) to be given by 

• 2^ ^ _ \ [TTW 1 

^"2 [3 1-/3^ 1 + ^{1 + 3P'){1-P')' 
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a > critical 




Fig. 5. 

Up to now we have considered a constant coupling a. We know that for q > X 
the Green function is determined by perturbation theory, which has to match 
the solutions in the region of smaller q. If > | for all g^, the solution which 

goes as ip ^ for q ^ matches the solution |('?/; — vr) ^ ^ + ^ for g — > oo 
monotonically. This determines mo as a function of rric in a unique way. Let A be 
the value of q where /3^(A^) = |. If, however, < | below q = X, the solutions 
can oscillate and mo(mcJ = for some m^. as indicated in Fig.6. 



m 




mo 



Fig.6 

This then is a solution corresponding to broken chiral symmetry. 
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